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Abstract 


We study the color and star formation rates of paired galaxies in filaments and sheets using the EAGLE 
simulations. We find that the major pairs with pair separation <50 kpc are bluer and more star-forming in 
filamentary environments compared to those hosted in sheet-like environments. This trend reverses beyond a pair 
separation of ~50 kpc. The interacting pairs with larger separations (750 kpc) in filaments are on average redder 
and low-star-forming compared to those embedded in sheets. The galaxies in filaments and sheets may have 
different stellar mass and cold gas mass distributions. Using a KS test, we find that for paired galaxies with pair 
separation <50 kpc, there are no significant differences in these properties in sheets and filaments. The filaments 
transport gas toward the cluster of galaxies. Some earlier studies find preferential alignment of galaxy pairs with 
the filament axis. Such alignment of galaxy pairs may lead to different gas accretion efficiency in galaxies residing 
in filaments and sheets. We propose that the enhancement of star formation rate at smaller pair separation in 
filaments is caused by the alignment of galaxy pairs. A recent study with SDSS data reports the same findings. The 
confirmation of these results by the EAGLE simulations suggests that the hydrodynamical simulations are 


powerful theoretical tools for studying galaxy formation and evolution in the cosmic web. 


Key words: methods: data analysis — methods: statistical — galaxies: evolution — galaxies: interactions — 


(cosmology:) large-scale structure of universe 


1. Introduction 


Galaxies are the fundamental units of the observed large- 
scale structures in the Universe. Understanding their formation 
and evolution is one of the primary goals of modern 
cosmology. The growth of primordial density perturbations 
via gravitational instability eventually leads to the formation of 
dark matter halos. Dark matter halos represent peaks in the 
density field. The halos are surrounded by a diffuse neutral 
hydrogen distribution after the recombination. They accrete the 
gas, a process which radiates away their kinetic energy and 
causes gas to settle down at their centers. The cooling and 
condensation of gas at the centers of these halos are believed to 
be the primary mechanism for galaxy formation (Silk 1977; 
White & Rees 1978). 

The formation and evolution of galaxies are expected to be 
influenced by both the initial conditions at the location of their 
formation and their interactions with the surrounding environ- 
ment. Galaxies interact with other galaxies in their neighbor- 
hood. The galaxy-galaxy interactions are known to amplify the 
star formation rate (SFR) in galaxies (Barton et al. 2000; 
Nikolic et al. 2004; Woods & Geller 2007; Ellison et al. 2010; 
Patton et al. 2011). The environments of galaxies have crucial 
roles in their evolution. The colors and SFRs of galaxies are 


strongly affected by the local density of their environment 
(Gómez et al. 2003; Lewis et al. 2002). The galaxies become 
redder and low-star-forming in the higher density environments 
(Kauffmann et al. 2004). The suppression of star formation can 
be driven by different physical mechanisms. Ram pressure 
stripping is a common phenomenon in galaxy clusters (Gunn & 
Gott 1972). Galaxies in high density regions are more likely to 
encounter harassment (Moore et al. 1996, 1998), starvation 
(Larson et al. 1980; Somerville & Primack 1999), strangulation 
(Gunn & Gott 1972; Balogh et al. 2000) and gas expulsion by 
supernovae, AGN or stellar winds (Cox et al. 2004; Murray 
et al. 2005; Springel et al. 2005). Star formation in a galaxy 
may also be quenched through several other routes. The mass 
(Bimboim & Dekel 2003; Dekel & Birnboim 2006), morph- 
ology (Martig et al. 2009), presence of a bar (Masters et al. 
2010) and high angular momentum (Peng 2020) can cause the 
star formation activity in galaxies to cease. 

Many other galaxy properties depend on the environment. 
Elliptical galaxies are more commonly observed in dense 
groups and clusters (Oemler 1974; Dressler 1980; Davis & 
Geller 1976; Guzzo et al. 1997; Goto et al. 2003). Spiral 
galaxies mostly occupy intermediate and low density regions of 
the Universe. These environmental dependencies of clustering 
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are reflected in different statistical measures such as the 
correlation function (Zehavi et al. 2005), genus (Park et al. 
2005), filamentarity (Pandey & Bharadwaj 2005, 2006), local 
dimension (Pandey & Sarkar 2020) and mutual information 
(Pandey & Sarkar 2017; Bhattacharjee et al. 2020; Sarkar & 
Pandey 2020). The environment of a galaxy is generally 
characterized by the local density. The local density undoubt- 
edly plays a decisive role in galaxy evolution. However, it 
cannot completely characterize the environment of a galaxy. 
Early generation redshift surveys revealed that galaxies are part 
of an all-inclusive network comprising clusters, filaments and 
sheets surrounded by vast empty regions (Joeveer & 
Einasto 1978; Gregory & Thompson 1978; Einasto et al. 
1980; Zeldovich & Shandarin 1982). Galaxies and their host 
halos are embedded in different environments of the cosmic 
web (Bond et al. 1996). Pandey & Bharadwaj (2008) find that 
star-forming blue galaxies trace a more filamentary distribution 
compared to red galaxies. More than 80% of the baryonic 
budget in the Universe is accounted for by low density gas 
(warm-hot intergalactic medium, WHIM) in filaments (Tuomi- 
nen et al. 2021; Galarraga-Espinosa et al. 2021). Consequently, 
the gas accretion efficiency of dark matter halos in different 
environments may differ in a significant manner. Thus, the 
cosmic web can have a significant impact on the galaxy 
properties and their evolution. The galaxies that are located in 
different parts of the cosmic web can experience different 
physical conditions, such as different densities of gas, different 
levels of tidal forces, and different frequencies of interactions 
and mergers. 

The interactions between galaxies with comparable masses 
are known as major interactions. Such interactions trigger new 
star formation in galaxies. Interacting pairs can be hosted in 
different morphological environments of the cosmic web. 
Galaxy pairs are more frequently observed in denser regions. 
The filaments and sheets, being the denser parts of the cosmic 
web, can host a significant number of major galaxy pairs. In a 
recent work, Das et al. (2023) analyze Sloan Digital Sky 
Survey (SDSS) data to compare the SFR and color of major 
pairs hosted in filaments and sheets. They find that the major 
galaxy pairs with separation «50 kpc are relatively high star- 
forming and bluer when hosted in filaments. Contrarily, the 
major pairs at separations larger than 50kpc show a 
significantly higher SFR and bluer color in sheet-like 
environments. This behavior may be related to the preferential 
alignment of galaxy pairs with the filament axis reported in a 
number of recent works (Tempel & Tamm 2015; Mesa et al. 
2018). Star formation in a galaxy is primarily regulated by its 
available gas mass. The inflows and outflows (Dekel et al. 
2009; Davé et al. 2012) can significantly modulate the gas mass 
in a galaxy. Transient events like interactions and mergers can 
drive the galaxies out of equilibrium. The alignment of galaxy 
pairs with filament spines may lead to anisotropic accretion and 
higher gas accretion efficiency in these galaxies. In this work, 
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we intend to verify these findings using hydrodynamical 
simulations. 

The Evolution and Assembly of GaLaxies and their 
Environments (EAGLE) simulation (McAlpine et al. 2016) is 
a hydrodynamical simulation that studies the galaxy formation 
and evolution in a cosmological volume. It describes the 
formation of galaxies by gas falling into dark matter halos and 
their subsequent cooling and condensation. It would be 
interesting to study the color and SFR in major pairs in 
filaments and sheets using EAGLE simulations. In observa- 
tions, galaxy pairs are usually identified by applying simulta- 
neous cuts on the projected separation and the velocity 
difference of the galaxies in the rest frame. However, all these 
pairs may not be undergoing interactions. Some of the pairs 
identified in observations may not be close in three-dimensions 
due to chance superpositions in high-density regions like 
groups and clusters (Alonso et al. 2004). Also, we cannot 
construct a mock catalog for the observational sample of galaxy 
pairs used in Das et al. (2023) due to the smaller volume of the 
EAGLE simulations. So, we decided to use the real-space 
positions of galaxies available in simulation to identify the 
major pairs. This would avoid any errors in identification of 
galaxy pairs due to the projection effects. We identify the 
geometric environments of galaxy pairs using the local 
dimension (Sarkar & Bharadwaj 2009). Our primary aim of 
this work is to study interaction induced star formation in 
filaments and sheets using EAGLE simulations. This would 
help us to assess the roles of filaments and sheets in galaxy 
evolution. 

We organize the structure of the paper as follows: we 
describe the data and method of analysis in Section 2 and 
present the results and conclusions in Section 3. 


2. Data and Method of Analysis 
2.1. EAGLE simulation data 


The EAGLE simulation (McAlpine et al. 2016) is a set of 
cosmological hydrodynamical simulations in periodic, cubic 
comoving volumes ranging from side of length 25 to 100 Mpc. 
Such a simulation tracks the evolution of both baryons and dark 
matter in the Universe from a redshift of 127 to 0. The 
simulation adopts a flat ACDM cosmology with Q, = 0.693, 
Qn = 0.307, NQ, = 0.04825 and Hy —67.77 km ! s ! Mpc"! 
(Planck Collaboration et al. 2014). 

We download the various properties of galaxies from the 
publicly available EAGLE run simulation. We extract informa- 
tion on the position of the center of mass of galaxies in three- 
dimensions within a comoving cubic volume of size 100 Mpc? 
from the Ref — L0100N1504. Subhalo table. We consider the last 
snapshot of the simulation having Snapnum — 28 which 
corresponds to redshift z — 0. We select only those galaxies 
which are flagged as Spurious = 0. This ensures that we select 
only the genuine simulated galaxies by discarding all the unusual 
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objects with anomalous stellar mass, metallicity or black hole 
mass. We also download the SFR, stellar mass and cold gas mass 
of simulated galaxies using the Ref — L0100N1504 Aperture 
table. These are estimated within a spherical three-dimensional 
(3D) aperture with radius 30 kpc centered at the location of the 
minimum gravitational potential of a galaxy. Use of this criterion 
gives well suited stellar mass and star formation estimates as 
compared to observational results and is also recommended for 
use by the EAGLE simulation team (McAlpine et al. 2016). We 
also consider only those galaxies with a stellar mass >0. 
Combining the two tables with GalaxyID, we obtain all of the 
above mentioned information for 325,358 galaxies. We also 
extract the rest frame broadband magnitudes of galaxies 
estimated in u and r band filters (Doi et al. 2010) from the 
Ref — L0100N1504. Magnitude table, where u and r respectively 
denote the Ultraviolet and Red filter bands of SDSS. We 
combine this table with Ref — L0100N1504. Subhalo and the 
Ref — L0100N1504 Aperture table using GalaxyID to get all 
the required information. The magnitudes of galaxies in different 
SDSS filters are also computed in 30 kpc spherical apertures 
(Trayford et al. 2015). Finally, we have all this information for 
29,754 galaxies. For the rest of the analysis, we refer to u—r color 
of galaxies as the difference in their rest frame non-dust 
attenuated absolute magnitudes in u and r band. Only the 
magnitudes of the galaxies with stellar mass log(Mgtetlar/Msun) > 
8.5 are provided in the Ref— L0100N1504 Subhalo table. 
However, here we use the stellar mass estimates of galaxies 
from the Ref — L0100N1504 Aperture table where the mini- 
mum stellar mass of a galaxy is log(Mstenar/Msun) ~ 8.2. 
Observations show that the galaxies with stellar mass 
Meotettar < 3 X 10 Msun are actively star-forming and the galaxies 
having stellar masses above this critical value are generally 
quiescent systems (Kauffmann et al. 2003). For the present 
analysis, we consider only those galaxies which have their stellar 
mass in between 8.5 < log(Mtettar/Moun) < 10.5. Our mass 
limited sample contains a total of 21,305 galaxies. 

We identify the nearest neighbor in three-dimensions for 
each galaxy in our mass limited sample. We denote the distance 
to the nearest neighbor for each galaxy by r. Here distance 
refers to the 3D physical separation between the center of mass 
of the galaxies. Initially, we label each galaxy and its nearest 
neighbor in our sample as a possible pair. We then select only 
those pairs for which r < 200 kpc. We also apply a cut on their 
stellar mass ratio 1 < = < 3 to include only the major pairs in 


our analysis. This provides us with a total of 2264 major pairs. 
The smallest pair separation in our sample is ~6 kpc. 

We determine the morphological environment of the 
galaxies in the EAGLE simulation by estimating their local 
dimension (Section 2.2). We use GalaxyID to cross match 
these galaxies with our pair sample. The cross-matching yields 
a total of 2537 galaxies in major pairs. We find that 373 and 
276 out of these galaxies are residing in filaments and sheets 
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Table 1 
This Table Shows the Definition of Different Geometric Environments Based 
on the Local Dimension (D) of the Galaxies 


Morphological environment Local dimension 


Di 0.75 x D « 1.25 
D2 1.75 < D < 2.25 
D3 D 22.75 

DLS 1225 € D < 1.75 


D2.5 2.25 SD « 2.15 


respectively. It may be noted that we cannot determine the local 
dimension of all the galaxies in the simulation. 


2.2. Geometry of the Local Environment 


We characterize the different geometric environments of the 
cosmic web using the local dimension (Sarkar & Bharad- 
waj 2009). The local dimension is a simple measure based on 
the galaxy counts within spheres of different radii centered on a 
galaxy. The galaxy counts within a sphere of radius R centered 
on a galaxy can be written as, 


N(«R) —A RP, (1) 


where D is the local dimension and A is an arbitrary constant. 
The number counts N( « R) would scale differently with the 
radius R depending on the local geometry of the embedding 
environment. The radius of the sphere around each galaxy is 
varied between R, Mpc < R < Rə Mpc and the galaxy counts 
are measured for each radius. Only the galaxies that have at 
least 10 neighboring galaxies within this range are considered 
for this analysis. We fit the observed number counts N( « R) to 
Equation (1) using a least-squares fitting. The goodness of each 
fit is determined by estimating the X? per degree of freedom. 
We only retain the fits having x < 0.5 (Sarkar & Pandey 2019) 
and discard the rest from our analysis. We choose Rı = 2 Mpc 
and R = 10 Mpc for the analysis presented here. 

The local dimension D describes the morphology of the 
embedding environment. Ideally, a filamentary environment 
should have D — 1 and a sheet-like environment should have 
D -—2. A homogeneous distribution in three-dimensions is 
represented by D — 3. However, the filaments, sheets, clusters 
and voids are not idealized structures and they can have a wide 
variety of shapes and sizes. Each geometric environment is 
assigned a range of local dimension as listed in Table 1. A 
nearly straight filament is represented by a D1 type environ- 
ment. Similarly, a D2 type environment represents a two- 
dimensional sheet-like structure. A D3 type environment is 
embedded in a 3D distribution with a homogeneous nature. The 
galaxies can also reside near the junction of different types of 
morphological environments. D1.5 and D2.5 can be treated as 
intermediate environments. 
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Figure 1. The left panel of this figure plots the cumulative mean u-r color as a function of pair separation for major pairs residing in D1 and D2 type environments. 
The right panel displays the cumulative mean SFR for the same pairs in two different environments. We use 10 Jackknife samples to estimate the 1c error bars shown 


at each data point. 


3. Results and Conclusions 


The cumulative mean of the u—r color for the major galaxy 
pairs as a function of the pair separation is plotted in the left 
panel of Figure 1. We compare the results for the major pairs in 
filaments and sheets in the same panel. 

This affirms that the major pairs with pair separation 
r<50kpe are on average bluer in a filamentary environment 
compared to those residing in a sheet-like environment. 
However, this trend only persists up to a pair separation of 
~50 kpc. A crossover of the two curves corresponding to D1 
and D2 type environments is observed at r ~ 50 kpc. The major 
pairs with pair separation r > 50 kpc are significantly redder in 
filaments compared to those located in sheets. We also analyze 
the SFR in major pairs residing in filaments and sheets and 
show the results in the right panel of Figure 1. We find that the 
major pairs at closer pair separation (<50 kpc) in filaments are 
comparatively more star-forming than those located in sheets. 
We see an exactly opposite trend for the major pairs with larger 
pair separation (>50 kpc). The colors of the galaxies are 
strongly correlated with their SFR (Baldry et al. 2009) and the 
results depicted in the two panels of Figure | are consistent 
with each other. It is also interesting that the crossover is 
observed at nearly the same pair separation (~50 kpc) for both 
color and SFR. We estimate the lo error bars at each pair 
separation using 10 Jackknife samples drawn from the original 
data sets. 

The stellar mass (Birnboim & Dekel 2003; Dekel & 
Birnboim 2006; Bamford et al. 2009) and the available cold 
gas mass content (Saintonge et al. 2012; Violino et al. 2018; 
Thorp et al. 2022) also play a very important role in deciding 


the SFR in galaxies. We test if the differences occurring in u—r 
color and SFR of galaxies in major pairs residing in D1 and D2 
type environments arise due to the differences in their stellar 
mass and cold gas content. We apply a Kolmogorov-Smirnov 
(KS) test to compare the distributions of stellar mass and cold 
gas mass of major paired galaxies in D1 and D2 type 
environments. The probability distribution functions of the 
two properties in D1 and D2 type environments are visualized 
in the two panels of Figure 2. We first carry out the test for the 
major pairs with all possible pair separations. We then conduct 
separate tests for the major pairs with pair separation >50 kpc 
and «50kpc. The results for the KS test are tabulated in 
Table 2. We note that the null hypothesis for all the major pairs 
can be rejected at the 90% and 99% confidence levels for stellar 
mass and cold gas mass respectively. This implies that the 
stellar mass distribution of galaxies in major pairs residing in 
D1 and D2 type environments is likely to be drawn from the 
same parent population. However, the galaxies in major pairs 
from filaments and sheets have a significantly different cold gas 
mass distribution. We also arrive at the same conclusions for 
the major pairs with r > 50 kpc. Interestingly, the results for the 
major pairs with r < 50 kpc suggest that the null hypothesis for 
stellar mass can be rejected at a very low confidence level 
(<60%), whereas for cold gas mass, it can be rejected at the 
<90% confidence level. Thus, stellar mass of major pair 
galaxies with r « 50 kpc in D1 and D2 type environments is 
highly likely to be drawn from the same underlying population. 
This clearly shows that stellar mass and available cold gas mass 
of the paired galaxies are not responsible for the differences 
observed in their u-r color and SFR in D1 and D2 type 
environments at smaller pair separations (r < 50 kpc). 
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Figure 2. The left panel displays the probability distribution function of log(Mgtettar/Msun) for major pairs in D1 and D2 type environments. The right panel features the 


probability distribution function of MaSScolagas for the same pairs. 


Table 2 
This Table Shows the Summary of the Kolmogorov-Smirnov Tests carried out for a Comparison of log(Mgtettar/Mgun) and MaSScolagas Of major pairs in filaments and 
sheets 
Dks Dxs(a) 
Major pairs log(Mgtetar/Meun) MasScoidgas(10° Mgun) 99% 90% 80% 7096 6096 
All 0.1037 0.5105 0.1292 0.0972 0.0852 0.0773 0.0712 
r « 50 kpc 0.1032 0.3512 0.3754 0.2826 0.2478 0.2249 0.2072 
r 2 50 kpc 0.1055 0.5104 0.1388 0.1043 0.0915 0.0830 0.0765 


Note. Separate comparisons are also carried out for major pairs having r < 50 kpc and r > 50 kpc. The table lists the KS statistic Dgs along with the critical values 
Dxgs(o) beyond which the null hypothesis can be rejected at a given confidence level. 


Each galaxy is believed to have formed within a dark matter 
halo. The properties of the galaxy are expected to be intimately 
connected to the mass of the dark matter halo. In fact, the mass 
of the dark matter halo is believed to be the most important 
parameter that determines the properties of a galaxy (Corray & 
Sheth 2002). The amount of substructures in dark matter halos 
increases with increasing halo mass (Gao et al. 2004; Pandey 
et al. 2013). There are observational evidences in favor of the 
correlations between substructure and star formation fraction in 
galaxy clusters (Bravo-Alfaro et al. 2009; Cohen et al. 2014). 
Substructures can also influence the stellar population in the 
galaxy (Helmi 2020). We depict the distributions of halo 
masses in the paired galaxies in filaments and sheets in 
Figure 3. The halo masses are obtained within the same 
aperture as the galaxies. We perform a KS test to find that the 
halo mass distributions of the paired galaxies in sheets and 
filaments are significantly different (Table 3). We ascertain that 
the halo masses of the paired galaxies in sheets are relatively 
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Figure 3. This figure displays the distributions of the halo mass of paired 
galaxies in sheets and filaments. 
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Table 3 
This Table Shows Kolmogorov-Smirnov Statistic Dks for Comparison of log(Mhaio/Msun) of Major Pairs Residing in D1 and D2 type Environments 


Dxs 
Major pairs log(Mnao/ Msn) 99% 
0.6532 0.1292 


Dxs(@) 
90% 80% 70% 60% 
0.0972 0.0852 0.0773 0.0712 


Note. It also shows the critical values Dgs(a) above which the null hypothesis can be rejected at different confidence levels. 


more massive than those residing in the filaments. The effects 
of the halo mass may also come from the virial shock heating 
of the halo gas that becomes important at masses greater than 
10^ Msun (Bimboim & Dekel 2003). Such heating can 
suppress the supply of cold gas by preventing cold streams 
from the intergalactic medium. However, we find that none of 
the paired galaxies in filaments and sheets in our sample reside 
in such a massive dark matter halo. At low masses, the 
supernova feedback may expel or heat the gas reservoir and 
quench the star formation (Kaviraj et al. 2007). The halo mass 
may have a role in shaping the physical properties of the galaxy 
pairs in filaments and sheets, but it is difficult to explain the 
crossovers observed in Figure 3 using these differences in the 
halo mass distributions. 

The filaments appear at the intersection of sheets and are 
generally denser compared to the sheets. Studies with N-body 
simulations suggest a successive flow of matter from voids to 
sheets, sheets to filaments and filaments to clusters (Ramachan- 
dra & Shandarin 2015; Galárraga-Espinosa et al. 2023). A 
number of earlier studies find that the galaxy pairs are 
preferentially aligned with the filament axis (Tempel & 
Tamm 2015; Mesa et al. 2018). The alignment signal is 
reported to be stronger for closer pairs residing near the 
filament spine. The anisotropic accretion along the filaments 
may significantly influence the gas accretion efficiency in these 
aligned galaxy pairs and trigger interaction. induced star 
formation in them. Contrary, the major pairs with 
r > 50 kpc show less star-formation in filaments than in sheets. 
The filaments are generally denser than the sheets. The D1 type 
galaxies are embedded in a high density environment as 
compared to the D2 type galaxies (Pandey & Sarkar 2020). The 
galaxies in denser environments are known to be redder and 
less star-forming (Lewis et al. 2002; Gómez et al. 2003; 
Kauffmann et al. 2004). So, naively one would expect the 
galaxies in a filamentary environment to be less star-forming 
and redder compared to the galaxies in a sheet-like environ- 
ment. We find that this is true for the galaxies in major pairs 
with separation larger than 50 kpc. However, the galaxies in 
major pairs at closer pair separation exhibit a strikingly 
opposite behavior. 

We do not analyze the alignment of the galaxy pairs in our 
study. The individual sheets and filaments cannot be identified 


using the local dimension. We plan to carry out a detailed study 
of the galaxy pair alignment with different identification 
techniques of the cosmic web in a future work. 

The EAGLE simulation provides two definitions for the 
position of galaxies. These are based on the center of mass and 
the location of the minimum of the gravitational potential. The 
two positions do not coincide for some galaxies. In this work, we 
use the center of mass to define the position of galaxies. We also 
repeat our analysis considering the minimum of the gravitational 
potential as the position of a galaxy. We show the results of this 
analysis in Figure 4. The main findings of our analysis remain 
unchanged with this alternative definition of galaxy position. 
Further, it is important to ensure that the major pairs considered 
in our analysis do not belong to the galaxy groups. We measure 
the distances to the 5th nearest neighbors for the paired galaxies 
in sheets and filaments and find that ~20% of them have their 
5th nearest neighbor within a distance of 500 kpc-1 Mpc. We 
discard these galaxy pairs and repeat our analysis. The results of 
this analysis are displayed in Figure 5. We find that discarding 
such galaxy pairs does not alter our results. 

The results reported in this work are very similar to the 
results obtained in a recent study (Das et al. 2023) of the color 
and SFR of major pairs in filaments and sheets using the SDSS 
data. Das et al. (2023) rely on a volume limited sample of 
galaxies (M, x — 19) for their analysis and find a crossover in 
these properties at nearly the same length scale (~50 kpc). It is 
interesting to note that we observe exactly the same trend in the 
EAGLE simulation data. This provides strong theoretical 
support to the observational findings that large-scale structures 
like sheets and filaments affect galaxy interactions. This also 
indicates that the galaxy properties are modulated by the 
geometry of their large-scale environment. 

Finally, we conclude that the filaments play a significant role 
in deciding the color and SFR in galaxies. The observed 
differences in the color and SFR of major pairs in filaments and 
sheets cannot be interpreted in terms of the differences in the 
local density and the stellar mass distributions. The interacting 
galaxy pairs with smaller pair separation can trigger star 
formation. The filaments provide a favorable environment for 
such interactions. This makes the interacting galaxies bluer in 
filaments compared to those found in sheets. 
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Figure 4. Same as Figure | but when the galaxy positions are defined by the location of the minimum of the gravitational potential instead of the center of mass. 
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Figure 5. Same as Figure 1 but after discarding the major pairs for which the 5th nearest neighbor lies within a distance of 500 kpc-1 Mpc. 
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